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Abstract 

An implementation of the multiple-scattering approach to x-ray magnetic 
circular dichroism (XMCD) in K-edge x-ray absorption spectroscopy is pre- 
sented. The convergence problems due to the cluster size and the relativistic 
corrections are solved using an expansion of the Dirac Green function for 
complex energies up to the second order in 1/c. The Fermi energy is dealt 
with via a complex plane integration. Numerical methods used to obtain the 
semi-relativistic Green function in the whole complex plane are explained. 
We present a calculation of the magnetic circular dichroism at the K-edge of 
bcc-iron including the core hole effect. A good agreement is found at high 
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energy. The physical origins of the XMCD spectrum near the edge and far 
from the edge are analyzed. The influence of the core hole, the possibility 
of a multiple-scattering expansion and the relation of XMCD with the spin- 
polarized density of states are discussed. A simple interpretation of XMCD 

at the K-edge is presented in terms of a rigid-band model. 
78.70Dm, 78.20Ls, 75.10Lp, 75.50Bb 
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I. INTRODUCTION 



X-ray magnetic circular dichroism (XMCD) in a magnetic sample is the difference be- 
tween absorption spectra obtained from right- and left-circularly polarized x-rays. This 
experimental technique was discovered in 1987,0 and gives information on the magnetic con- 
tribution of each orbital and of each atomic species in a sample. Still now, the detailed 
mechanisms that give rise to XMCD at the K-edge are not entirely clear near the edge and 
are unknown far from it. The purpose of the present paper is to understand XMCD at the 
K-edge over the whole energy range used in experiments. Our main tool is a new semi- 
relativistic equation which is very useful from a computational point of view, and which can 
be used to calculate other relativistic properties, such as magnetic anisotropy. 

In a previous paper, we described a preliminary approach to the calculation of the x-ray 
magnetic circular dichroism effect in x-ray absorption spectroscopy, within the framework of 
the multiple-scattering theory.l We assumed that XMCD at the Ln.m-edges of rare earths 
and transition metals has a simple interpretation and that XMCD at the K-edge was the 
most difficult case, because spin-orbit acts on the photoelectron, which makes a physical 
picture difficult to build. /^From the experience that was gained in the intervening years,il 
we know now that the Ln in-edges of rare earths are more complicated than expected, and 
some progress has been made towards a reliable practical use of the experimental resultsi'i 
at the K-edge of transition metals. 

Although fully relativistic calculations of XMCD at the K-edge are numerous,iHii we 
believe that the use of a semi-relativistic approach is still justified. Its main advantage is that 
spin and space variables are uncoupled at the zeroth order. This enables us to use the full 
local cluster symmetry, to build a clearer physical picture, and to make orientation averages. 
The smaller scattering-matrix dimension of the semi-relativistic approach allows for the 
calculation of bigger clusters and broader energy ranges. Finally, the existing multiple- 
scattering numerical programs can be easily extended to calculate the XMCD. For instance, 
Ankudinov and Rehr have adapted their FEFF program^ to calculate the extended structure 
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in XMCD at the Ln,iii-edge of rare earths and 5d transition metals. They obtained good 
agreement for gadohnium. A similar adaptation is possible, although more complex, for 
K-edge spectra. 

Here we report on the calculation of the K-absorption spectra and the XMCD of Fe in 
the presence and absence of a static core hole in the Is state. The influence of the core hole 
was taken into account by the final-state rule which assumes that the final state energies 
of the x-ray absorption process are measured in the presence of a static core hole, i.e., the 
dynamics of the excitation process are neglected.iil This rule has been used with success 
to explain the satellite structures on the high energy side of the emission spectra of simple 
metals.illii In particular, for Na and Al the Ln ni emission satellites on the high-energy side 
of the main lines due to the double ionization of the 2p core level is well reproduced by 
the calculations if the excited atom is treated as an impurity and the static electron-hole 
interaction is included in the calculation.illii The final state rule was also used to calculate 
the Lin-edge of 3d transition metal ferromagnetics,^ and the effect of the core hole is found 
to be very important. The determination of the self-consistent potentials is carried out 
within the local density approximation with a linear muffin-tin orbital (LMTO) basis-set.ii 
We use a supercell geometry and treat the excited atom with a core hole as a single impurity 
atom in a lattice. A Is electron is added to the conduction states in the supercell, and we let 
the system of N-l-1 electrons relax selfconsistently. We have found that the photo-electron 
localizes on the excited atom providing an efficient screening of the core hole. Like the 
case of Ln,iii edgeil up to 80% of the photoelectron polarizes as a minority spin. Thus, the 
electronic structure of the excited Fe atom is very close to that of a cobalt metal. However, in 
this calculation we have found that the effect of the core hole on the K-absorption edge and 
the XMCD is small, and this is because the p-states of Fe are less affected by the presence 
of the core hole. 

In this paper, the determination of the absorption spectra and the XMCD is carried out 
within the multiple scattering theory. During the implementation of our early formalism we 
encountered several difficulties that convinced us that the naive approach used in Ref. ^ to 
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take spin-orbit effect into account is matfiematically not sound. In tfie present paper, we use 
a different matfiematical framework for the calculation of XMCD, that can be generalized 
to take relativistic effects in the valence d-shell for Lnjn-edges. This approach gives a 
direct expansion of the Dirac Green function in powers of ^, it overcomes the divergence of 
the Foldy-Wouthuysen transformation through the use of complex energies, and it leads to 
results in reasonable agreement with experiment. 

Some of our results are useful only for x-ray absorption spectroscopy, but others can have 
wider applications, such as the convergence of the semi-relativistic limit and the properties 
of the Green function at complex energies. 

The remainder of the paper is organized as follows. In the second section we point 
out the limitation of the Foldy-Wouthuysen transformation and give a new method for the 
determination of the relativistic effects within the Green function approach. We then apply 
the Green function framework to the calculation of the XMCD signal. In a third section the 
use of complex energies is discussed, and a complex plane integration is used to take the 
Fermi energy into account. In the fourth section we present the numerical implementation 
that we used, and give the new XMCD cross-section formula at the K-edge. Finally, we 
present a calculation that was carried out on a reasonably large Fe cluster of 259 atoms 
near the edge and a cluster of 51 atoms up to 500 eV. Good agreement with experiment is 
found in the high-energy region. Various issues, such as the influence of the relativistic core 
hole, the origin of XMCD at the edge and far from it, the possibility of a multiple-scattering 
expansion are discussed. The spectra arc related to the density of states through a simple 
rule, and the rigid-band model is found to be correct at high energy. An appendix showing 
how symmetrized bases were used to calculate XMCD faster and to understand some aspects 
of the experimental spectra closes the paper. 



II. THE RELATIVISTIC EXPANSION 



5 



A. The Foldy- Wouthuysen transformation 



The relativistic corrections to physical phenomena are usually treated using a perturba- 
tion approach developed in 1950 by Foldy and Wouthuysen (FW),il who made a unitary 
transformation to eliminate the coupling of the large small components from the Dirac equa- 
tion and obtained, to zeroth order, the Pauli Hamiltonian and, to first order in (1/c^), the 
spin-orbit, Darwin and kinetic energy corrections. Later, this Hamiltonian was taken as a 
starting point for the calculation of many physical properties, such as magnetic anisotropy, 
XMCD or magneto-optic effects. Sometimes, the spin-orbit term is used in a second-order 
perturbation theory, overlooking that the FW-transformation must then also be made up to 
second order (1/c^) to achieve formal consistency. However, when carrying out higher-order 
transformations, strongly divergent terms are found even for the simple case of a hydrogen 
atom:@ the second order "correction" is a sum of infinite terms. Additional information 
concerning the convergence properties of the FW-transformation can be found in Ref. p7 . 

Even at the 1/c^ order, the behavior of the Hamiltonian obtained by adding the kinetic 
energy correction to the non-relativistic Hamiltonian is so bad that all negative eigenvalues 
(bound states) disappear, and a continuum of states is obtained between — oo and mc^/4, 
however small c may be.ii In some cases,ilii only the spin-orbit interaction is added to the 
non-relativistic Hamiltonian. For a Coulomb potential, the negative spin-orbit term pulls 
the particle into the singularity and the resulting Hamiltonian is not essentially self-adjoint 
(it can have any eigenvalue, depending on the chosen self-adjoint extension) Physically 
reasonable results could be obtained by choosing the boundary condition where the wave 
function is zero at r = but this process is mathematically ambiguous.il 

To summarize, the series obtained from the FW-transformation exhibits three kinds of 
divergence: (i) the series itself diverges, which is usually not a problem since most pertur- 
bation series used in quantum mechanics are divergent (e.g. Zeeman effect. Stark effect, 
anharmonic oscillator) ;!! (ii) from the second order term (1/c^), each term of the series is a 
sum of infinite terms that should add to a finite term, but no procedure is known to carry 
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out this summation;E3 (iii) the first order term of the series (spin-orbit couphng) diverges as 
the cluster size increases. 

Recently, all these problems of the semi-relativistic limit were solved using modern math- 
ematical methods, and we propose to use this new formulation for the calculation of XMCD. 
Other approaches to the semi-relativistic limit of the Dirac equation were proposed by quan- 
tum chemists.S'ii 



B. Relativistic (1/c) expansion of the Green function 



The standard method to calculate relativistic effects within the Green function approach 
is to consider the spin-orbit and other terms as a perturbation of the Hamiltonian, and to use 
the Lippmann-Schwinger equation corresponding to that perturbation.iil This procedure 
is not mathematically safe and leads to divergences. Another approach to XMCD was 
proposed by NatoliU where the Schrodinger equation including spin-orbit was solved for 
each muffin-tin sphere. 

Here, we start from the retarded Green function G^{z) corresponding to the Dirac Hamil- 
tonian with potential V (the Dirac Green function) and we use a slight modification of the 
analytic 1/c expansion of G^{z) that was obtained in Ref. ^ (see Ref. ^ for a detailed 



proof): 



-1 



G{z) G{z)Ay{2mc) 
AG{z)/{2mc) [AG{z)A^ + 2m]/{2mcf 



(1) 



where 



Tiz) 



/ 

V 



/ 



(2) 



G{z)A^{V - z)/{2mc) 
[AG{z)A^ + 2m]{V - z)/{2mcf 

All matrix entries are 2x2 matrices. A = —iha.V and a are the Pauli matrices and G{z) 
is the Green function for the Pauli Hamiltonian with potential V. V is generally a 2x2 
matrix that describes the potential experienced by the electrons with up and down spins. 
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In the simplest case, V is a diagonal matrix made up of and V^. The successive terms of 
the relativistic expansion are obtained by the series [1 — T{z)]~^ = 1 + T{z) + T^(z) + ■ ■ -. 
This expansion gives a very compact formulation of the correction terms. For instance, the 
first order correction to a non-degenerate bound-state eigenvalue Eq corresponding to the 
eigenstate \ip) is given by 

AE = {^\A^{V-Eo)A\^), 

H which contains the kinetic energy correction, the Darwin correction and the spin-orbit 
correction. Because of this simplicity, formal manipulation can be carried out much further, 
for instance to obtain a consistent second order expansion of the spin-orbit interaction. 
Another aspect of this expansion is that the small components are not eliminated, which 
will be useful for x-ray absorption spectroscopy. 

It was shown that expansion (|l]) is analytic,^ but its actual radius of convergence was 
investigated only recently. Preliminary results00 for the Dirac Green function without 
potential showed that the radius of convergence was a function of the imaginary part of 
the energy. For a Dirac equation with potential, the Green function expansion could be 
shown@ to be convergent only when the imaginary part of the energy is greater than mc^ (a 
very bad resolution indeed). However, the convergence can be much better when a physical 
property is calculated instead of the general Green function. For instance, the convergence 
of the bound state energy expansion is fast Physically, the idea that can be drawn from 
these mathematical results is that the use of complex energies makes the semi-relativistic 
expansion converge, whereas it diverges on the real axis. This idea will be detailed in the 
next section. 

/^From the formal point of view fully relativistic programs require a coupling of an infinity 
of different / values due to the magnetic field.ii Up to now, the coupling between / and / ± 2 
was neglected, which limits the maximum possible value of /.@ In our approach, this coupling 
is made consistently with 1/c through the expansion process. 
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C. Application to XMCD 



In x-ray absorption spectroscopy, the spin-orbit parameter for the core hole has a mag- 
nitude of several hundreds of eV and cannot be considered a priori small. Furthermore, rel- 
ativistic effects in x-ray absorption have been shown to be considerable both theoretically^ 
and experimentally.^ Therefore, a fully relativistic core hole wave function must be used. 
This is not possible within the standard semi-relativistic approach. In the electric dipole 
approximation, the relativistic formula for the x-ray absorption cross-section is: 

a(e) = -47rao/icu(i|(e*.r)Im[G'^(r,r';2)](e.r')|i), (3) 

where z = huj + Ei + ie, Ei is the energy of the core level \i) and e is the x-ray polarisation 
vector. The calculation of the fully relativistic initial state \i) presents no difficulty, thanks 
to the availability of Desclaux' program.0 In Ref. ^ we supposed that the relativistic effects 
on the valence and continuum states are weak. We show in the present section that this 
supposition is not always true. 

The influence of relativistic effects on the photoelectron is obtained by the series expan- 
sion of the Dirac Green function. The Dirac Green function in Eq.(|I]) is expanded up to 
second order in 1/cS 

{i\e*.vG^t.Y'\i) ~ (0|e*.rGe.r'|0) 

+ 7^{<P\e*.vGA^e.v'\iP) + ^{^\e* .vAGLv'\<p) 
2mc 2mc 

+ j^^{<l>\e*.rGA^{V - z)AGe.r'\<jy) 

+ j^^^{tP\i*.r{2m + AGA^)iy\^) (4) 

where the two-component spinors |0) and {ip) are the large and small components of the 
Dirac spinor \i). In our previous approach! we have neglected the relativistic nature of the 
core state, and used only the spin-orbit interaction, as a consequence we had obtained only 
the first and a part of the fourth term. Formula is valid for a non-muffin-tin potential; 
the use of non- muffin-tin potentials is probably required to achieve a quantitative estimate 
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of XMCD. In all the terms of Eq. (^) only the fourth term yields magnetic circular dichroism 
for a powder sample. Therefore, in the following, we call ctq the spectrum obtained from the 
first term, 

ao(e) = -AnaohujY^{<t>^\ie\r)lm[Giv,r';z)]ie.r'M^) (5) 

s 

where the sum over s is the sum over initial states (two spin states for a K-edge). For a 
powder, ctq does not depend on e. We call ai the spectrum obtained from the fourth term 
of Eq. (I). 



cri(e) = -A7Taohuj'^{(j)''\{i* .r) Im 



{e.v'm (6) 



and <7mcd = )~cri(e+) for an external magnetic field aligned with the x-ray wavevector. 
More precisely, we have in the fourth term:S 

A\V - z)A = -n'^liVV.V) + {V- z)A + ta.{VV x V)]. (7) 

In the right-hand side of the above equation, only the third term contributes to XMCD at 
the K-edge, because the first two terms do not connect space and spin variables. In the case 
of spherical potentials, one obtains the usual formula for the spin-orbit interaction:^ 

We took this formula for convenience, to make contact with Ref. |], but it is also possible 
to act with A and A'' directly on the Green function. This would be much simpler if all 
relativistic effects (and not only XMCD) were investigated.^^ In principle, o"i(e~) +cri{e^) 
contributes also to the normal absorption spectrum, but this relativistic contribution was 
assumed to be small as compared to ctq. 

Finally, to obtain the influence of all sites, the potential V must be written as a sum 
over all sites plus an interstitial potential. Because the operator A is translation invariant, 
each site can be taken as the origin when making transformation (^. In the muffin-tin case, 
it can be shown that the interstitial region does not contribute to XMCD. 
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For real energies, XMCD calculated within the present approach is formally identical 
with the results of Ref. 0. Therefore, we can use the latter to understand the real energy 
behavior of the former. For small clusters, the calculation of XMCD was easy,0 but as 
the cluster size increased, we met difficulties. Fig. ^ shows ctq and aucD for a cluster of 
259 atoms of bcc-Fe. Both spectra exhibit a very sharp resonance (width 0.005 eV). 
This resonance is due to the fact that, as the cluster size increases, one comes close to the 
singularities of the crystal scattering matrix.0 The second point is that (Tmcd is larger than 
ctq- This unphysical result proves that the spin-orbit effect is not a small perturbation of the 
system in that energy range. This illustrates the fact, discussed in the previous section, that 
the relativistic expansion is generally not convergent at real energies. It may be worthwhile 
noticing that this divergence of the relativistic expansion is not due to the Z/ singularity 
of the spin-orbit interaction in a Coulomb potential. 

Moreover, the experimental resolution at the K-edge of transition metals is of the order 
of 1 eV whereas Fig. ^ shows much sharper structure. Therefore, the usual method of con- 
voluting the theoretical spectrum with a Lorentzian profile for comparison with experiment 
would mean calculating much more points than actually needed. 

For real potentials, the singularities of the multiple-scattering matrix are on the real axis, 
and the use of a complex energy smoothes the sharp structure of (Tq linked with the cluster 
size. We show in the next section that the convergence problem of the relativistic expansion 
is also solved by calculating the Green function at complex energies. 

Finally, we work within a real space multiple-scattering approach, so that the cluster is 
assumed to be of finite size. The relativistic expansion of the Green function for an infinite 



crystal was studied in Ref. 52 



III. COMPLEX ENERGIES 

In this section, we justify the use of complex energies, and we show how this simplifies 
the calculation of physical quantities. We assume that the potential is a diagonal 2x2 
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matrix made up of and V^, which are the potentials experienced by up and down spins. 
These potentials are assumed to be real. The generalization to complex potentials is a non- 
trivial task,ii and the use of an energy- dependent width T(E), which can be considered to 
be the imaginary part of the potential (constant over space), was found to be sufficient to 
reproduce experimental spectra.ii 



A. Green function 

The link between the wave function and the Green function formulae for x-ray absorption 
cross-section is established through the identityiil 

E \f){mE - Ef) = 6{E -H) = -^[G^E) - G-{E)], (9) 

where G^(-E') = G{E ± ie), e is a positive number that is taken to tend to zero at the end 
of the calculation (this is just a notation to designate an integration path in the complex 
plane) . 

Because of the finite core-hole lifetime, the calculated spectrum must be convoluted by a 
Lorentzian with a half width at half maximum (HWHM) F (which may ultimately depend 
on the photoelectron energy). In the Green function formalism, the convolution with a 
Lorentzian is obtained by calculating the Green function for a complex energy.!! This can 
be shown from the formula: 

r+oo TG^(e)-G'(e) 

L (E- eP + r' - + - - ('o' 

Since G(r, r'; z) = G{r', r; and G{r, r'; z)* = G{r', r; z*) for Hermitian potentials V, we 
have G{E - iV) = G{E + iV)* . 

Moreover, the term huj = e — Ei oi the absorption cross-section can be taken care of 
by noticing that zG{z) = HG{z) + 1. Thus eG'^^e) = HG^{e) + 1, which gives, after 
convolution, HG{E + iT) + l = {E + iT)G{E + iV). 

Therefore 



12 



ao(e) = -47r«o(«|(e*.r) lm[{E + zT - E,)G{r,r'; E + iT)]{i.r')\t). 



(11) 



In other words, the convoluted spectrum is obtained by calculating the Green function 
for an energy with an imaginary part F. 



B. Fermi energy 



In the previous section, we did not take into account that all one-electron states up to 



the Fermi energy Ep are occupied. To yield physical results, integral (|T0D must be carried 
out from Ep instead of from — oo. This modifies our expression in an interesting way. We 
must now evaluate : 



-/ 

2i Je 



de 

1 r+oo 

4:71 jr.p 



V{e-E{){G+{e)-G-{e)) 



{E - e)2 + F2 

/ de{e-E,){G+{e)-G~{e)) 



(12) 



.E-e + iV E-e-iTi 

This integral can be calculated by a complex plane integration technique:!^ Since the 
self-consistent potential V is Hermitian, the poles of ^^(e) are at e = r — ie, where r is 
a real number. Therefore, we can choose the contour of Fig. |^ to apply Cauchy's integral 
formula.ii The contribution of the Jordan contour at infinity is zeroii and we obtain: 



]_ T {e-E,)G+{e) 
2iJep ""tt (E-e)2 + F2 



--e{E - Ef){E + tT- Ei)G{E + tT) 



F j-oo (Ep + tt - Ei)G{EF + It) 
\jo {Ep + it- Ey + 



(13) 



Heavyside's step function 6{E — Ep) is present in the expression because, when E > Ep, 
the pole E + iT is inside the contour. 

The second integral in the rhs of Eq. (0) is closed on the negative imaginary side, and 
we obtain finally: 

F(e-Ei)Im[G+(e)] 



de 



F 

— / dtRe 

TT JO 



vr {E- e)2 + F2 

\Ep + it - Ei)G{Ep + it) 
{Ep + it- E)2 + F2 



e{E - Ep)lm[{E + iT- Ei)G{E + iT)] 



(14) 
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It can be shown that the right-hand side of Eq. ([M]) is continuous at E = Ep in spite 
of the step function.0 The recourse to this contour integration is efficient, as compared 
to an integration over real energies, because the Green function is quite smooth on the 
fine z = Ep + it (see Fig. ^) and only a few points must be calculated to evaluate the 
integral. Moreover, G{Ep + it) tends rapidly to the local atomic Green function of the 
absorbing site, as shown in Fig. |. Therefore, G{Ef + it) can be obtained at high t just by 
calculating a simple one-site Green function. The vanishing of the neighbor's influence can 
be understood through the form of the structure constant matrix H (see section The 
if -matrix elements have the form H"^^ = exp {iK,Rij)P{l/ uRij) where P is a polynomial. For 
large t, k = y/z ^ ^t/2{l + i), and the if-matrix elements describing the influence of the 
neighbors are damped by a factor of exp —^t/2Rij. The integration in the complex plane 
was chosen along a straight line. This is not a steepest descent path but was found to be 
sufficiently efficient. 

The above proofs are the same whether G{z) is the non-relativistic or the relativistic 
Green function (assuming that the negative energy states are full). To obtain the expression 
for o"! after convolution, we start from the results obtained with the Dirac Green function 
and use the relativistic expansion (^). This gives an expression for ai which is the sum of 
Eq. (P) for z = E + iT and of an integral over the line z = Ep + it. The detailed form of 



the result will be given in section |1VC . 

It was observed by Rehi§ that thermal motion can also be a convergence factor in 
multiple-scattering calculations. Although thermal effects can be formally accounted for 
in full multiple-scattering calculations,^ temperature effects were neglected in the present 
study. 

Finally, a characteristic of the 1/c^ relativistic term is that spin- flip is not allowed for 
powders. An explicit calculation of the spin-flip scattering amplitude shows that it is an 
order of magnitude smaller than non-spin-flip scattering for Nd^"*".^ Experiments^ conflrm 
that elastically scattered electrons are rarely spin-flipped. 
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IV. NUMERICAL ASPECTS 



The x-ray absorption cross-section is obtained from the cluster Green function for com- 
plex energies, which we write asHS 

1^ —Ti Wml'm' + -0t/'0m,m'0i,j) (15 

where k = ^/z, 6\ is the (complex) phase-shift for potential \^*(r), t\ = siia6}expi6l, R\{r) 
is the regular solution of the radial Schrodinger equation for potential \/*(r) that matches 
smoothly to cos5^j^(/tr) — sin5^n^(/tr) at the mufhn-tin radius pj, Hl{r) is the irregular solu- 
tion of the radial Schrodinger equation for potential \^*(r) that matches smoothly to h^{nr) 
at the muffin-tin radius. Finally the multiple-scattering matrix is r = [T~^ — kH] ^, where 

{TaYL,^, = -{ty^)kA^'Sm,m' and HZ,^, = -47rzEA,^'-^"-''Qi™A>^('^i?.,)>T(%)- In 
the last expression, the Hankel function is defined as the function h^^ of Ref. |6^. This 
straightforward extension of the real energy case is possible because wave functions and 
phase-shifts are analytical functions of 

We present here the numerical methods that were used to calculate the Green function 
and the ground state. First, we describe how the radial wave functions were obtained, then 
we show how the Green function singularities were avoided, and finally we explain how the 
multiple-scattering matrix was calculated. 



A. Radial wave functions 

The core level energy and wavefunction were obtained from the Dirac SCF program of 
Desclaux.0 

The photoelectron wave functions were obtained by adapting the method of Ref. ^ to 
complex energies. A non-normalized regular wave function is written as ^(r) = r^fe{r) 
with fe{r) = fr{r) + ifi{r). For the complex energy E = Ej. + iEi the radial Schrodinger 
equation is : 
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= {V{r) - Er)fr 

= (\/(r) - E^)U 



r 

2(f + 1) 
r 



f'i - -£^i/r, 



(16) 



with the boundary conditions /^(O) = 1, /;(0) = -Z/{1 + 1), /^(O) = 0, //(O) = 0. Where 
/' and /" are, respectively, the first and second radial derivatives of /. This system of 
second order differential equations is transformed into a system of four first order differential 
equations which is solved by a forward step-adaptative fourth-order Runge-Kutta method.il 
The phaseshifts 5^ are deduced from the Wronskian equation : 



exp {2i5£) 



(17) 



hj{Kp)g^{p) - Khj {kp)R^{p) 

The normalized radial wave functions R\{r) are obtained from gi^ir) and the phaseshifts. 

The irregular wave function is written as = r^^^^^^g^^r) with gi{r) = gr{r) + igi{r) 
and obtained by backward step-adaptative fourth-order Runge-Kutta method^ from the 
boundary conditions iff(p) = hf{Kp), H/{p) = Khl\np). 



B. Multiple-scattering matrix and singularities in the Green function 

Formula (62) of Ref. ^ could not be directly applied because of possible infinite normal- 
ization factors arising from the matching of the wave function at the muffin-tin radius to 
a sum of Bessel functions.0 With the present normalization at the muffin-tin radius, the 
radial wave functions cannot be singular, but sin 5^ terms arise in the denominator of the 
Green function (see Eq. p!5|) . Therefore, we must examine the behavior of the scattering 
matrix when sin 5] becomes small.S We have 



r = TJI - KHTaY' =Ta + KTa[l - KHTa]'^ HT^. 



(18) 



Hence, 



^ij^^/'f^w + ([1 - i^HTa] ^H)^ sin(5;^,e 



(19) 



16 



and we see that r|^^„/ sin^ 5^ is generally singular at sin^^ = because of the first term. 
These singularities did not appear in previous calculations of x-ray absorption cross section 
because Im(r) was used, which is regular since an additional sin 5^ factor comes from e**^^. 
To avoid these singularities, one can use (T^mfm + sin^ (J^, which is smooth. Therefore, 

all terms of the Green function ( ]TB| ) are regular. 

To avoid the singularities discussed in the previous section, and to use a more symmetric 
equation, we use a polar decomposition technique^ and define a modified multiple scattering 
matrix f by r = + \fT'af\fT'a. Numerically, f is obtained from the equation : 

f= [l-z^V^if^T;]"'-!. (20) 

f is a regularized multiple-scattering matrix that describes the effect, on each atom, 
of the rest of the cluster (f is zero for a cluster of one atom). The matrix elements of 
[1 — K^T'aH\/Ta\ are calculated from the phaseshifts and from the efficient algorithm for 
the calculation of H that was proposed by Rehr and Albers0 and tested against alterna- 
tive methods.^! The form [1 — K\fT'aB. \/T'c\ ^ combines the advantages of the two standard 
equations r = [T^^^ — kE.\ ^ and r = Ta[l — kHTo^^ . As in the first equation, a sym- 
metric matrix is inverted (when real spherical harmonics are used), so that fast inversion 
algorithms for symmetric matrices can be used.0 As in the second equation, it is regular 
when iTa)^lrniimi ~ 0' ^'^ ^'^^^ spurious singularities of are very simply avoided.^! More- 
over, asymptotic analysisH shows that it is well behaved when i is large and/or k is small, 
whereas none of the two standard equations is. The matrix is inverted using a standard LU 
decomposition technique.!! 

C. Cross-sections 

The detailed expression of the absorption and XMCD cross-sections for a powder relies 
on the orientational averaging technique used in Ref. 0, and more calculation steps are given 
in Ref. As in Eq. (|I^) , all cross sections are written as the sum of a Green function term 
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and an integral. With our modified multiple-scattering matrix, the "Green function" part 
of the K-edge absorption cross sections at energy E = hu + Ei with a HWHM T becomes: 

ao = E + zr) + a^jE + tV)] (21) 

s 

where the sum over s is the sum over the two spin states (i.e. the sum of the cross- 
sections calculated for potentials and V^). For each spin state one defines aoa{z) = 
{AnaQ/ 3) {z — Ei)iy/z exp[i6i{z)]D^ (z) , which is the atomic contribution to x-ray absorption, 
and 

aoniz) = {47rao/3){z-E,)V~zD'{z)exp[zS'^{z)] ^ . ' ' ^ (22) 

v3 smdi[z) 

which describes the influence of the neighbors {i = denotes the absorbing site). We recall 
that! f^^iU, aa; z) = (-l)(i-^')(l-/i£m|aa)ff^,„(2;). 

The radial integrals are D{z) = /q°° r^dr(f)o{r)R^{r; z), where 0o(^) is the large component 
of the core hole wave function and D^{z) = r^dr(f)Q{r)F{r; z) where F{r, z) is an auxiliary 
function defined as F{r;z) = r'^dr' (j)Q{r')R\{r^] z)H^{ry] z) {R\ and HI are defined in 
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The "Green function" part of the magnetic circular dichroism cross-section is written: 
aMCD = Im[E((-l)(-V2)(^.^(^ + ,r) + ai,{E + tT) + al^{E + zT))] (23) 

s 

where aiai^z) = (47rao/3)(-2 — Ei)z exp[2i6l{z)]M^^ (z) describes the purely atomic contri- 
bution to XMCD (the Fano effect), 

an{z) = -(47rao/3)(z - E,)2tz exp[2t5',iz)]D{z)M'' (z) ^ , \ ] / (24) 

v3smdj'(z) 

is the local contribution due to the spin-polarization of the p-states on the absorbing 
site, and 

ai„(z) = (47rao/3)(^ - E,)zD\z) ^ exp[z(<5?(^) + 5l{z))]a{z) 

X E [{i - a){i + a + 1) + 2] E(-l)°'" "°'^^^:""' (25) 
a=\i-i\ a sm6^{z)smdi{z) 



describes the contribution to XMCD due to the spin-orbit scattering of the pho- 
toelectron by the neighbors and the absorber. The radial matrix elements are 
M^{z) = S^Pdri{r)Rl{r-z)F{r;z) , M""{z) = j;^ r'^dr^{r)F\r; z) and Ci{z) = 

2 

r'^dr^^ {r)[Rl{r; z)] , with ^{r) as defined in Eq. (^). 



The analysis of § [iV B| shows that f^^{li, aa; z) can be divided by y sin 5i{z) sin 5l{z)^ and 
all terms of the XMCD cross-section are now regular. 

For each term, the presence of the Fermi energy is taken into account by calculating 
integrals along the line e = Ep + it as in Eq. (|14|) : to all \m\a{E + iV)] terms we add the 
integral: 



- / dtRe 

TT JO 



a{EF + it) 



{Ef + it- Ef + r2 



(26) 



V. RESULTS 

The converged potentials used in the multiple scattering formalism were obtained by 
means of an all-electron self-consistent, scalar-relativistic and spin-polarized linear-muffin- 
tin orbital (LMTO) method.ii We used the exchange-correlation potential and energy in the 
von Barth and Hedin approximation. For the Brillouin zone integration of the density of 
states we used the tetrahedron method with about 300 k points in the irreducible part of 
the Brillouin zone.0 To simulate the effect of the core hole, we treated the excited atom as 
a single impurity in a lattice using a supercell geometry. We have used increasingly larger 
supercells to ensure the convergence of the magnetic moment and the density of states of 
the impurity site. The final calculations were done for a simple cubic lattice of 16 atoms per 
unit cell; the lattice parameter being a = 2ao- The results for the magnetic moment and 
the density of states are close to those of the supercell of 4 atoms. To use larger clusters 
in the multiple-scattering formalism we had to assume that the potentials for distant shells 
are bulk like. 
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A. Near edge region 



Fig. and b show the results we obtained for a converged cluster of 259 atoms of bcc-iron 
(diameter 2.0 nm). The convergence was investigated by checking that the spectral shape 
becomes stable with respect to the cluster diameter and by comparing to the results obtained 
with a cluster of 821 atoms (diameter 2.9 nm) on a wider energy mesh. The potentials for 
the initial state (without core hole) and for the final state (with core hole) were obtained 
by a self-consistent super-cell calculation as indicated above. We used touching muffin-tin 
spheres without overlap. Because of the core hole width (F = 2 eV), absorption is possible 
at energies lower than the Fermi level. The effect of the Fermi level is clear, especially for 
the XMCD spectrum. With a Fermi level Ep = 2 eV, the first positive structure disappears 
whereas with Ep = —2 eV, it is too large. Therefore, the size of the first positive peak 
depends strongly on the position of the Fermi level. 

The absorption and XMCD cross-sections are now smooth and the divergences due to the 
cluster size and the relativistic corrections have disappeared. We see also that the spectra 
are continuous at the Fermi energy, although they are given by two different formulas for 
E > Ep and E < Ep. This corroborates the fact that our treatment of the Fermi level is 
numerically sound. 

Figure ^ compares our calculation with experiment for the edge region. All experimental 
and theoretical spectra were normalized so that the absorption edge jump is 1. The normal 
spectrum is not well reproduced. The agreement for the XMCD spectrum is better, although 
the excellent agreement in the intensity of the first two peaks is fortuitous, because the degree 
of circular polarization was not one for the experimental spectrum. For the calculated 
spectrum, various contributions are presented. The full thick line is the total contribution, 
and the full thin line is the local contribution {an of Eq. (p^). In Fig. ain as written in 
Eq.(p5D is expanded into the contribution of each i (£=1,2,3) and j=0,l,2 (absorbing atom, 
first shell, second shell). For later use, the results are given without taking Fermi energy 
into account. No term i=0 exists because £=0 gives no spin-orbit. For the absorbing site. 
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the terms with even i are zero because of symmetry. /^From Fig. it can be observed that 
the first positive peak (peaking at 1 eV) is not due to the local term (0), but to the spin- 
polarization of the (i-states of the neighbors (mainly the term corresponding to 1=2 and j=l 
in Eq. (|25D). This is related to the weak ferromagnetic nature of bcc-Fe,@ and was already 
observed on a smaller clusteiS and in tight-binding calculations.0@ Both absorption and 
XMCD calculated structures are too large from 20 eV above the edge. This reduction, 
which is very common in multiple-scattering calculations, is probably due to photoelectrons 
that experience inelastic interaction with the metal. This inelastic effect comes into play 
above the plasmon energy (about 10 eV) and can be included in our calculation through an 
energy-dependent The main failure of the XMCD calculation is the presence of a large 

second positive peak which is absent in the experiment. A similar peak can be observed in 
fuUy-relativisticil and tight-binding calculations, although a direct comparison is difficult 
because of the different normalization used. 



B. Core hole 

Two modifications of the core hole were tested. In the first one, the non relativistic 
(Schrodinger) equation was solved for the core hole. The core hole energy was found different, 
but the Is wavefunction was quite similar to the relativistic one, and the normal spectrum 
was a bit smaller (because the normalization of the relativistic wavefunction includes the 
small component, which does not contribute to the spectrum) but the XMCD spectrum 
was not distinguishable from the relativistic result. Similarly, core hole exchange splitting is 
negligible since the difference between the core states obtained with up and down potentials 
did not yield noticeable effects. 

The second test was conducted to test the influence of the core hole on the spectra. 
Fig. 1^ shows the normal and XMCD spectra with and without core hole and with a Fermi 
energy Ef=^ eV). The normal spectra are quite similar, except for an overall amplitude 
factor. The XMCD spectra are more different, although the position of the structure does 
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not move. Comparison with experiment does not enable us do decide for a model. 

C. SPEXAFS 

In XMCD at the K-edge, spin-orbit acts directly on the photoelectron, and it is interesting 
to know whether the high-energy part of the spectrum is related to the spin polarization of 
the p-states of the absorbing site. Using the methods developed in Ref. ^ it is possible 
to show that, in the single scattering approximation, only the p-projected term o"|; — aj; 
survives at high energy. Figure H(b) shows the spectrum obtained with a cluster of 51 
atoms, calculating all angular momenta up to £=8, and multiplied by 1/3. The cluster is 
too small to be realistic in the edge region, but above 30 eV, the overall agreement is correct. 
The peak at 110 eV in the theoretical spectrum should be broadened, and the peak at 60 eV 
in the experimental spectrum is probably due to multielectronic effects, which are known to 
be strong in XMCD spectra at that energy.H The thin line represents the contribution of the 
p-projected states, which is seen to be dominant at high energy. Moreover, the calculated 
XMCD spectrum is in phase with the calculated EXAFS spectrum. This phase relation 
between dichroic and normal spectra was observed experimentally by Pizzini and coll.B 

It is sometimes assumed that XMCD reflects the spin-polarization of p states projected 
on the absorbing atom (p^ — p^). Fig. |^ shows that the absorption spectrum is indeed quite 
similar to the p density of states, but the XMCD spectrum is quite different from the spin 
polarization of the p density of states. In fact, the spin-polarization is very similar to the 
derivative of the density of states. In other words, the rigid band model becomes correct at 
high energy (a somewhat surprising result), and the band splitting is about 1 eV. 

/^From this result, a very simple approximate expression can be derived for XMCD. 
Assuming a rigid band model, we can consider that the up and down bands are exchange- 
splitted by the energy AE. Moreover, neglecting the non-diagonal terms of the spin-orbit 
operator, one can consider that the m = ±1 components of the p-band are split by spin-orbit 
coupling (i.s. Therefore, for transitions towards i = l,m = 1 final states (left-circularly 
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polarized x-rays), we have 

a+T = a{E + - Q (27) 

= a{E - AE/2 + C) (28) 
CT+ ~ 2a{E) + - CfSa/dE'^ (29) 

and for transitions towards £ = l,m = — 1 final states (right-circularly polarized x-rays): 

a-^ = a{E + AE/2 + () (30) 
cr-i = a{E - AE/2 - () (31) 
a- c:^2a{E) + {AE/2 + Cfd^cr/dE^ (32) 

therefore, XMCD becomes 

a+ - a" ~ -2AECd^a/dE\ (33) 

Since C is fairly constant at high energy, the exchange-splitting energy AE can be deduced 
from experimental spectra. Figure |TU| shows XMCD together with the second derivative of 
the normal spectrum (multiplied by 10). Because the model is very crude, the agreement 
is not perfect, but good enough to say that the image deduced from the rigid-band model 
is correct. The presence of this second derivative explains also the phase relation between 
EXAFS and XMCD (considering EXAFS as a sum of sines). A more elaborate interpretation 
comes from considering the expression for ao„ (Eq. ([2^)) and an (Eq. (p^). Both expressions 
involve the multiple-scattering matrix f°°(ll, 00; 2), but an contains an additional factor i 
that makes it proportional to the real part of f^^{ll, 00; z), whereas the EXAFS cross section 
is proportional to its imaginary part (it was checked that the other factors do not intervene 
much in the phase at high energy). Since real and imaginary part of the Green function are 
related by Kramers- Kronig theorem, which transforms sine functions into cosine functions, 
and because dichroic effect is due to the difference between spin up and spin down an, 
it is proportional to the derivative of the real part of 'f°°(ll, 00; z) which, because of an 
additional minus sign, is in phase with the imaginary part of the multiple-scattering matrix 
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(to see this, consider that, in the EXAFS regime, the Green function can be approximate by 
exp(i/ci? + 0), the imaginary part is a sine function - the EXAFS formula - and the derivative 
of the real part is a sine function as well). Of course, this correspondence between XMCD 
and EXAFS is not exact, and much interesting information comes from the difference. 



D. Multiple-scattering expansion 

The program that we used to calculate the XMCD spectrum of iron is relatively heavy 
and slow. To make the analysis of XMCD spectra as easy as that of normal absorption 
spectra, it is necessary to investigate the validity of the multiple-scattering expansion which 
is used by much faster programs, such as FEFF.S To do this, we compare the results obtained 
by the full inversion of (1 — KTaH)~^Ta with the single scattering expansion (T^ + nTaHTa + 
K^TaHTaHTa) and the double-scattering expansion (previous term plus K?{TaH)^Ta). This 



comparison is shown in Fig. |ri|. The overall agreement is correct but not excellent, probably 
because of the shadowing effect which is large in bcc structures. Notice that, in principle, 
the zeroth-order scattering (T^ + nTaHTa) also contributes due to the term cri„. However, 



this contribution is very small at high energy. Higher order terms are considered in Ref. ISB. 



VI. CONCLUSION 

Since the first experimental XMCD spectra,0 a number of calculations of XMCD K-edge 
spectra have been carried out, using different methods. A fully relativistic KKR method 
was used for bcc-Fe,&l00 hcp-Co,lili Fe-Co alloys, Slfflllilli and cfc-Ni.0 A relativistic 
LMTO-ASA calculation of the K-edge XMCD spectrum of Fe in GdFe2 was carried out 
in Ref. 0, a molecular orbital approach was used for Fe in tetrahedral and octahedral 
environments,^ a multiplet approach for Ni in a molecular magnet@ and a tight-binding 
method for metallic Fe and NiB@ and Co,0 where XMCD was related to the projection of 
the orbital momentum along the x-ray direction. All these calculations were restricted to a 
narrow energy range around the edge. We have developed in this paper a multiple-scattering 
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approach which allows for the calculation of extended structure, and includes the core hole 
without additional effort. 

We have presented a solution of the convergence difficulties associated to the one-electron 
calculation of physical properties related to spin-orbit interaction. A recourse to Green func- 
tions with complex energy argument led us to smooth absorption and XMCD cross- sect ions. 
The presence of the Fermi level was accounted for through a complex plane integration 
which was found much more stable than on the real line. This technique can be used to 
calculate other spin-orbit influenced properties, such as anisotropy energy or spin- dependent 
spectroscopies. 

Robust and accurate numerical methods were proposed to evaluate the Green function 
in the whole complex plane, and the smooth behavior of the cluster Green function for large 
imaginary energies was explained. 

We saw that the EXAFS part of XMCD at the K-edge is simply connected to the spin- 
polarization of the p-states on the absorbing site. The high-energy part is therefore simpler 
to interpret than the near-edge part, where spin-orbit interactions with the neighbors gives 
strong contributions. 

Our first application is encouraging, and further comparison with experiment, including 
the relation between magnetic and non-magnetic fine structure,0 will be presented in a 
forthcoming publication. Multiple-scattering calculation of XMCD at the Ln,iii-edges of Gd 
are presented in Refs. ^ and ^7\. 

As observed by Ankudinov and Rehr,0 XMCD in the x-ray range provides a very good 
test of effective potentials representing the exchange interaction of the photoelectron with 
matter. An alternative exchange potential was proposed by Zhogov et al.0 in their treatment 
of magnetic EXAFS. In the present study, we have used the potential provided by the ground 
and relaxed-excited states, but did not include any specific exchange potential. 
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APPENDIX: SYMMETRIZED BASIS 

One of the advantages of the semi-relativistic limit of the Dirac Green function is that 
one is allowed to use the full local point symmetry of the absorbing atom without taking spin 
direction into account, since the spin and space functions are not coupled. In this section, 
we show how symmetrized bases were implemented to reduce the size of the problem. 

For each spin state, the cluster potential is written as V{r) = Z]j^'('^i) (* runs over 
atomic sites), and let Pa represent the action of the symmetry operation a on the cluster. 
If a belongs to the local point symmetry group Q of the absorbing atom {i = 0), then 
Pa[V(r)] = V{r). Therefore, the Green function G{z) = {z — Hq — V)^^ is also invariant for 
operations of Q. This can be used to reduce the size of the multiple-scattering matrices.ii0 

In practice, we set up a symmetrized basis by the following procedure. For an irreducible 



26 



representation (irrep) a of the symmetry group, we take a matrix realization V'^°'^ of a, and 
we define the (pseudo-)projector: 

y a&g 

where the cluster symmetry operator P^{a) has matrix elements in the \nim) basis : 

[P^{a)]n'm'nm = -D^/^ (a)5n',a(n) • (A2) 

Here, g is the number of elements of the symmetry group Q, da is the dimension of irrep a, 
P'm'mi^) is the (eventually improper) Wigner rotation matrix corresponding to the symmetry 
operation a of the point symmetry group We choose a column jo of the matrix realization 
(jo = 1) and, for each \nim) representing a spherical harmonics Yf^ attached to site n, 
we calculate the projection Pj^j^\nim) . Let p be the set of atoms that can be obtained 
from a given atom by operations of Q. For each i and each p one can generate a vector 
space E{ap£jo) by calculating all the projections of |n£m), for n G p and m = —i, 
By singular value decomposition,!! an orthonormal basis of E^apijo) is obtained, which is 
denoted by |ap£jos), where s = 1, dim{E{ap£jo)). This symmetrized basis can be written 
as a function of the initial basis kets : 

\apijos) = ^ {n£m\apijos)\nim) . (A3) 

m=-i,...,e 

The partners of \apijos) are obtained by 

\apijs) = P^flapijos) (A4) 

The symmetrized basis is then used to simplify the matrix inversion yielding the relevant 
multiple scattering matrix elements t°*^,. It can be shown that 

{P^{a)H)nlmn'l'm' = {H P^ {a'))nlmn'l'm' (A5) 

Therefore, 



27 



{apejs\H\a'p'e'j's') = {aplj,s\PlfHP\?:^f\a'p'ej',s') (A6) 

= {o^P^3os\Pi:]'P^^jH\a'p'l'3,s') (A7) 

= Sa,a'Sj,j'{apejos\H\ap'e'jos'). (A8) 

Similarly, the atomic scattering matrix is diagonal: 



where = sm6^ exp{id^) and is the £-th phaseshift of the atoms of type p. Because 
of the transformation properties of matrices H and T^, the full multiple-scattering matrix 
T = [T~^ — kH]~^ is diagonal in a and j: 

{apejs\r\a'p'£'j's') = Sa,a'Sjj'r^i]p,i,,,. (AlO) 

In other words, the multiple-scattering matrix can be inverted separately for each irrep and 
it repeats itself dim(r("^)) times in each irrep. For the example of a cluster of 259 Fe atoms, 
only the Ti^ irrep is relevant, because of dipole selection rules, and the matrix to be inverted 
has dimension 288 instead of 4144. Once r^lpnigi is obtained, the matrix elements which are 
required for the calculation of XMCD at the K-edge are obtained through the basis change: 

ajss' 

where pi is the class of atoms to which atom i belongs. 

An additional advantage of the use of a symmetrized basis is the fact that the inverse of 
r is more rapidly calculated. In practice, the matrix elements {ap£i s\H\ap' £' j s') must be 
calculated from the matrix elements H'^^,^, and the basis change matrix by 

{ap£jQs\H\ap'ejos')= ^ E («io^^^«l^^^)^^m^'m'(^'^'^'l«io^ys')- (A12) 

"^^P mm' 

Instead of summing over all pairs n e p and n' e p\ we can choose a member uq of class p, 
and a set of members n\ of class p' which span all possible non-equivalent neighbors of no 
in class p' . In other words, for each pair (n, n'), there is an element a of symmetry group Q 
and a representative rii such that a{n) — Uq and a{n') — rii. Therefore 
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We have also 

Y^DlM{nim\c^Jo(^Ps) = Y.^l,{a-'){noif^\ajips), (A14) 

m j 

which yields, using the fact that the symmetrized matrix elements of H do not depend on 
Jo- 

{apijos\H\ap'i'jos') = ;^ E E rj„,(a)rj„^.,*(a) (A15) 
X E(«J^Ps|noM^Sw(^i^'"^'l«j'^Vs')- (A16) 

mm' 

/^From the orthogonality theorem of group representations,!^ the sum over j" does not 
depend on a, and 

\n\N 

{apejos\H\ap'i'jos') = E ^-^^^-^ {aj ips\n,^^^ (A17) 

where \p\ is the number of elements of class p, and Nngm is the number of elements n' of 
such that a(no) = no and a{ni) = n', where a & Q. Therefore, the sum is reduced to pairs 
of inequivalent neighbors. 

/^From these symmetry considerations, one can show that, at the K-edge of a transition 
metal in a cubic environment, the spin-polarization of the (i-shell of the absorbing atom 
cannot be measured, but the spin-polarization of the (i-shell of the neighbors can. This 



was observed in the multiple-scattering approach in Refs. ^ and ^ and in the tight-binding 
approach in Ref. [f8|,[f^. 
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FIGURES 

FIG. 1. Calculated K-edge absorption spectrum (thick line) and x-ray magnetic circular 
dichroism (thin line) of a cluster of 259 atoms of bcc-Fe near the singularity at the vicinity of 
Fermi level. 

FIG. 2. Pole structure of the integrand in the second term of Eq. (p^). 

FIG. 3. (a) : Y^s^Oa + ^On ^^'^ (b) : Y.si~^Y~^^'^^ta + '^'i; + (^tn foi" complex energies 
along the line Ep + it {Ep = 0). In each case, the corresponding atomic quantities J2s^0a ™d 
J2s{~^y^^^'^^ta ^1'^ plotted, to show the convergence of the full Green function to the atomic 
Green function. 

FIG. 4. Calculated (a) K-edge absorption spectrum and (b) x-ray magnetic circular 
dichroism (XMCD) of a cluster of 259 atoms of bcc-Fe with different values of Fermi level, 
Ep = — oo, —2, 0, 2 eV. For the XMCD signal, the size of the first positive structure depends 
strongly on the position of Fermi level, and disappears for Ep above 2 eV. 

FIG. 5. Comparison of the calculated (a) K-edge absorption spectrum and (b) the x-ray 
magnetic circular dichroism (XMCD) of a cluster of 259 atoms of bcc-Fe with the experimental 
results of Ref. ^ The Fermi level is Ep=0 eV. 

FIG. 6. Contribution to ain of the absorbing site and the first two shells (from top to bottom) 
and for £=1,2,3. No Fermi energy was used. The contribution of a given shell is the contribution 
of one atom of the shell multiplied by the number of atoms in the shell. 

FIG. 7. Calculated (a) K-edge absorption spectrum and (b) x-ray magnetic circular dichroism 
(XMCD) of a cluster of 259 atoms of bcc-Fe with (thick line) and without (thin line) core hole. 
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FIG. 8. Comparaison of the calculated (a) K-edge absorption spectrum and (b) x-ray magnetic 
circular dichroism (XMCD) of a cluster of 51 atoms of bcc-Fe with the experimental results of Ref. ^. 
The Fermi level is Ef=0 eV, the core-hole broadening is F = 2 eV, and the maximum scattering 
wave is ^max=8. 

FIG. 9. (a) Comparison of the calculated K-edge absorption spectrum with the density of p 
states projected on the absorbing atom and (b) comparison of the calculated XMCD with the spin 
polarization of the density of p states projected on the absorbing atom. The derivative of the DOS 
is also presented, to illustrate the rigid band picture that becomes valid at high energy. 

FIG. 10. Comparison of the XMCD spectrum with the second derivative of the absorption 
spectrum (multiplied by 10) for a cluster of 51 iron atoms. 

FIG. 11. (a) Comparison of the calculated K-edge absorption spectrum with the single and 
double scattering expansion (b) comparison of the calculated XMCD with the single and double 
scattering expansion. 
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